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We briefly review some examples of confinement wliich arise in condensed matter physics. 
We focus on two instructive cases: the off-critical Ising model in a magnetic field, and an 
array of weakly coupled (extended) Hubbard chains in the Wigner crystal phase. In the 
appropriate regime, the elementary excitations in these 1+1 and quasi-one-dimensional 
systems are confined into 'mesons'. Although the models are generically non-intograblo, 
quantum mechanics and form factor techniques yield valuable information. 
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1. Introduction 

The phenomenon of quark confinement hes at the center of modern particle 
physics. What is perhaps less widely known is that similar phenomena also 
emerge in condensed matter physics. In this brief review, we discuss some 
examples of 'quark' confinement in 1+1 and quasi-one-dimensional systems. 
Celebrated cases include the Ising model in a magnetic field [1] and the 
spin-1/2 Heisenberg chain with frustration and dimerization [2,3]. As is 
frequently the case in one-dimension, Lorentz invariance may emerge as an 
approximate symmetry of the low energy effective action. This opens the 
way to an assortment of field theory methods [4] and to a wealth of analytic 
results. 

The 'quarks' in these, and many other cases, appear as domain walls 
(or solitons) interpolating between different vacua. In particular, these ex- 
citations may carry 'fractional' quantum numbers which differ from those 
of the electron. A prominent example are the excitations of the half-filled 
Hubbard model, where 'spinous' carry spin-1/2 and no charge, and 'holons' 
carry charge ibe and no spin. The existence of these fractional excitations has 
been confirmed by numerous experiments on materials composed of weakly 
coupled chains, e.g. [5,6]. 

Whilst these exotic excitations may be stabilized through the judicious 
choice of materials and parameter regime, they may form bound states in 
the presence of additional perturbations. In particular, the omnipresent 
interchain interactions may confine the elementary excitations and radically 
alter their masses and quantum numbers. A simple of illustration of this idea 
arises in an array of weakly coupled quantum Ising chains. The Hamiltonian 
governing a single chain is given by: 

H\\ = -j\\Y.{^i^Ui+9^f), (1) 

i 

where > is a ferromagnetic exchange constant, and the operators a^'^ 
are Pauli matrices residing at site i. As is well known, this model undergoes 
a r = quantum phase transition at the point g = 1 [7]. In particular, for 
g < I, the ground state degeneracy is spontaneously broken and the system 
develops long-range order {(t? ) ^ 0. Let us now consider the effect of a weak 
interaction between neighboring chains: 

Here, J± <C Jii , and for notational simplicity we consider a two-dimensional 
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array of chains. In the ordered regime, the effects of neighboring chains may 
be treated in a mean field approximation: 

^r = -E^<.; h=h^J^{an, (3) 

where Z± is the transverse coordination number of the lattice. In this man- 
ner one obtains decoupled Ising chains in an effective magnetic field: 

^ = E [-"^ll + 9 - h at] . (4) 

i 

We emphasize that this is not an external magnetic field, which one may 
of course apply, but it arises quite naturally through the interchain interac- 
tions. As follows from the pioneering work of McCoy and Wu [1], this weak 
magnetic field acts as a linearly confining potential on the zero field excita- 
tions: the 'quarks' are confined into a rich spectrum of 'mesons'. Moreover, 
in this regime of weak confinement, analytic progress is possible. 

The layout of this contribution is as follows. In ^11 we discuss the Ising 
model in a magnetic field in a little more detail. We focus on the properties 
in the ordered state where bound states form. In ^ we discuss an exam- 
ple of a quasi-one-dimensional Wigner crystal: a weakly coupled array of 
quarter-filled extended Hubbard models. This system has much in common 
with both the Ising [1] and Heisenberg [2] examples. In the limit of weak 
confinement, many 'generations' of 'mesons' are shown to exist. We conclude 
in SI 

2. Ising Model in a Magnetic Field 

The Ising model has had a venerable history and there are excellent reviews 
devoted to it — see for example [8,9]. In particular, the confining aspects 
of the 2D Ising model in a weak magnetic field were first exposed by McCoy 
and Wu [1]. The implications for gauge invariant correlation functions in 
the 2D Z2-gauge Z2-IIiggs system were subsequently investigated [10]. For 
a brief overview of these ideas see also page 106 of ref. [11]. 

In the ensuing discussion we contrast the behavior of the Ising model 
in the absence and presence of a magnetic field. In both cases we cast 
our discussion in terms of (suitably defined) form factors. The relevance 
of form factor techniques in this and other non-integrable models was first 
emphasized by Delfino, Mussardo and Simonetti [12]. 
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2.1. T T- and H = 



As is well known, in the scaling region close to criticality, the 2D Ising 
model is described by the field theory of a free Majorana fermion — see for 
example [13]: 

^FF = ^ y [^^* + + i^^^] ■ (5) 

The fermion mass m measures the departure from criticality, and we take 
m > in the 'low-temperature' phase we are interested in. The energy and 
momenta of these particles are conveniently parameterized in terms of the 
rapidity 9 

V = msinh^, £ = mcosh.9, (6) 

where h = c = 1. In particular, they form an highly efficient basis in which 
to compute correlation functions [14-17]: 

(0|c7(r,2:)o-(0,0)|0) = 

E / fl'{n , (0|^(r,x)|gi...g.)(g....gi|a(0,0)|0).(7) 

_m J-oo V^^j 



n=0,2,. 



Here (r, x) are Euclidean coordinates and a{T,x) is the continuum version 
of the lattice spin. In the ordered regime, the spin field only couples to 
intermediate states with an even number of particles.'^ Equivalently, 

( I cr(r, x) o-(0, 0) I ) = 

where 

n n 

P = m sinh 6i, E = m cosh 6i. (9) 



i=l i=l 



The matrix elements of o"(0, 0) (or any other local field) between the vacuum 
and the multiparticle states, are known as form factors. The computation of 
form factors is central to all integrable models [14-17] and the lowest order 
contributions yield valuable information about the long distance correlations. 



'^Tlie corresponding disorder operator /x{r, x) couples to an odd number of particles. By Kramers- 
Wannier duality the situation is reversed in the disordered 'high-temperature' phase. 
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In particular, the zero-particle form factor of the spin field is the spontaneous 
magnetization: 

(0|cj(0,0)|0) = (a) =?ni/8s, s = 2^/^^e-^/^A^/\ (10) 

where A = 1.28243 ... is Glaisher's constant. Likewise the two-particle form 
factor is well known [14]: 



( I a(0, 0) I ^1^2 ) = i ( ) tanh y ^ ) ' ^^'^^ 

The higher particle form factors are also known [14] but they need not con- 
cern us here.^ Substituting 1)10(1 and l|ll() into the expansion ((HJ one obtains 

(0|cj(r,x)f7(0,0)|0) = 
where 9± = {9i ziz 92)/2. Performing the integral over 9^, 

|/_.\|2 roo 

(0|cj(r,x)cj(0,0)|0) = |(cj)|2 + ^^ / d9 tanh^ 9 Ko(2rch 9), (13) 

where f = mr and r = Va;^ + r^. The integral in (|13|) may be evaluated in 
closed form. This yields the famous result of Wu and collaborators [18]: 

{0\a{T,x)a{0,0)\0) =\{a)\'g{r), (14) 

where the scaling function 

^(^) = 1 + ^ - Kg(f)] - rKo{f)Ki{r) + ^K2(f)^ + . . . (15) 

bears a remarkable connection with the solution of the Painleve III equation 
[18,19]. The large distance asymptotics of the spin-spin correlation function 
are readily extracted: 

(0|a(T,x)a(0,0)|0) ~|(c7)|2|l + i-^ + O(e-^0}; r^oo. (16) 

This example helps convey the efficiency of the form factor approach. 

In condensed matter applications it is useful to introduce the so-called 
dynamical susceptibility 

X{^^, k) = Xe{^, k)c:,^s-ioj, (17) 



''As in all integrable models, this involves solving the Form Factor Axioms [17] (or generalized 
Watson equations) with the appropriate S-matrix; in this case S = — 1. 
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where XE{^,k) is nothing but the Fourier transform of the Euchdean spin- 
spin correlation function 

/oo roo 
dx / dre^'^^-*'=^(T^cj(r,x)a(0,0)). (18) 
-oo J —oo 

Tr denotes time ordering. The definition H17() includes the analytic contin- 
uation to real frequencies and e is a positive infinitesimal. The dynamical 
structure factor, as measured by inelastic neutron scattering, is extracted 
from this: 

S{uj,k) = -lmx{uJ,k). (19) 

In this way, the two-particle contribution to the dynamical structure factor 
may be obtained 



where w > 0. We plot this in Fig. ^ 



0.1 



3 
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Figure 1 

of a magnetic field: 



Dynamical structure factor of tiie Ising model in the ordered regime and in the absence 
VIIP~W and A = |(<j>|2/m^- 



As is well known, the dynamical structure factor is a direct reflection 
of the elementary excitations. Indeed, the absence of any sharp peaks in 
Fig. ^ reveals that they are deconfined - the spin operator couples to pairs 
of excitations and a two-particle continuum exists above threshold. 
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In the next section we recall the dramatic change of this picture in the 
presence of a weak magnetic field [1]. 

2.2. T ^ T- and H > 

c ^ 

In the presence of a magnetic field the scaling region may be described by 
the non-integrable'^ Ising Field Theory: 

AwT = Aff + h j (fxa{x), (21) 

where ^ff is the action of free massive fermions. As follows from McCoy 
and Wu [1], the effect of the magnetic field in (|21|) is to confine the massive 
free fermions or 'quarks' into bound states or 'mesons'. In the limit of small 
magnetic fields, the masses of these bound states follow from the Schrodinger 
equation of two particles of mass m subject to a linear confining potential: 

-m^^ + ^'^'^(^)=^^^''^^)' 
where we employ the reduced mass. In this example, the confining parameter 
A or 'string tension' is related to the magnetic field by [1,22] 

X = 2h{a), (23) 

and the corresponding bound states have masses given by 

mb = 2m + £b. (24) 

As follows from (|2H) A has dimensions of [mass]^ {h = c = 1) and equation 
((221) is consistent on dimensional grounds. In accordance with the more 
recent work of Fonseca and Zamolodchikov [22] , equations (|22() and l|23p may 
be derived from the action of H21|) on an appropriate two-particle bound 
state. Their method is rather general and employs the known finite size 
form factors of the perturbing field (cr) in the unperturbed (free fermion) 
model. In particular it permits a systematic study of the relativistic Bethe- 
Salpeter corrections to equation and the resulting bound state mass 
spectrum. That is to say, the Schrodinger equation ((221) is a small momentum 
approximation which holds for small magnetic fields or weak confinement. 
For our present discussion however, equations (|22|) and (|23jl are sufficient. 

The Schrodinger equation (|22() is easily solved [23]. In the region x > 
we introduce the change of variables 

^ = {m\)^/^{x-£/\) (25) 



"^At T = Tc the model is actually integrablo and yields the celebrated Eg mass spectrum [20, 21]. 
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SO as to yield the Airy equation 



^^(0 = 0. (26) 



The normahzable solutions may be written in terms of an Airy function [24]: 

^{0 « Ai(0. (27) 

As expected, this wavefunction is oscillatory in the classically allowed region 
X < £/\ and damped in the classically forbidden region x > £/X. Since 
the potential A|x| is an even function, the wavefunctions have a definite 
symmetry and must be matched smoothly at the origin. In particular, the 
antisymmetric wavefunctions must vanish at x = 0. This leads to the energy 
level quantization condition: 



Ai 



m \ 1/3 



0. (28) 



The bound states of the Ising model in a weak magnetic field are therefore 
indexed by the zeroes of the Airy function [1]. Their masses are given by 



m 



= 2m + — Zfe; Ai(-Zb) = 0. (29) 



Equivalently, in notations closer to those of McCoy and Wu [1] 

mfe = 2m+^^^^^^^^^; Ji (U;,) + J_ i (Uf,) = 0, (30) 

where we have used ()23|) and Zi, = (Ah/2)^/'^. The absence of the symmetric 
wavefunctions from the known Ising model mass spectrum ()3U() is apparent, 
and we shall return to this point in the next section. As A — > there 
is a proliferation in the number of bound states and their masses densely 
populate the interval between 2m and 4m; in general, when the mass of a 
bound state exceeds twice the mass of the lightest 'meson' 2mi they become 
unstable. As A ^ we confine our attention to bound states in the vicinity 
of threshold 2m. 

Having discussed the spectrum of the Ising model in a weak magnetic 
field, we now turn our attention to the spin-spin correlation function. The 
confining magnetic field yields a spectrum of bound states of mass m;,. On 
general grounds we expect the spin-spin correlation function to have the 
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form: 

{0\a{T,x)a{0,0)\0)x 



, 1 J -oo ^'^ 



where |\I'6(^)) is an asymptotic state (to be discussed in the next section) 
describing the non-trivial bound state. For a general operator O of spin-s, 
we expect 

(0|O|^6(^)) = (0|O|^fe)e"^ (32) 

Since o" is a spinless operator, the matrix elements p2() are independent of 
rapidity. Performing the integral over rapidity yields: 



rib 

{0\a{T,x)a{0,0)\0)x = | ( a )|2 + - V |( | ct(0, 0) | Ko(mbr); (33) 

vr ^-^ 

6=1 

where Uh is the number of stable bound states. The main problem of course 
is that we do not know the rapidity independent matrix elements appearing 
in ()33() . let alone their multiparticle extensions. After all, the Ising field 
theory (|2H) is non-integrable, and the powerful axiomatic approach to in- 
tegrable models [17] does not apply. It is however a perturbation of a well 
understood integrable model. In the next section we shall illustrate how 
such single particle bound state matrix elements may be calculated directly 
in the limit of weak confinement. Understanding their detailed form is an 
obvious task in a more systematic approach to non-integrable models with 
confinement. Before embarking on this pursuit however, let us proceed a 
little with the benefit of hindsight. In particular, let us focus on the long 
distance asymptotics of ()33|) . We know that as the confinement parameter 
A — > we must recover the asymptotics H15() of the unperturbed system. 
Indeed, this important cross-check was performed in the seminal work of 
McCoy and Wu [1]. In addition, let us assume that as A ^ the matrix 
elements are the same for all bound states; we shall justify this below. In 
this manner, the long distance asymptotics of the connected contribution to 
read: 

( I air, X) aiO, 0) | )^,^o ~ E (^4) 



b=l 

As A ^ 0, the bound state masses ()29|) become closer together. Put another 
way, the bound state poles in momentum space must coalesce to form a 
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branch cut [1]. In this limit one may therefore convert the sum in H34|) to an 
integrah 

/ —;= ~^ / "0 exp 



y/rubr y/2mr Jo 



2^/m 



(35) 



where we have used the fact that the zeroes of the Airy function have the 
hmiting form [24] 

w [37r/2(6- 1/4)]^/^ 6 ^ oo. (36) 
This yields the long distance behavior 

(0|.(.,.).(0,0)|0)5^.-.KMI^£^^ (3r) 

In order to recover the correct asymptotics (jlSp as the confining perturbation 
is turned off, we require that the bound state matrix elements have the 
specific form: 

|(0|a(0,0)|*,)|=i^./|; (38) 

Equivalently, 

|(0|a(0,0)|1'b)| = (o-)\/^; h ^ 0. (39) 

where h = {a)h/m?. Indeed, substituting ()39() into H34() we recover the 
expected asymptotics of McCoy and Wu [1]: 

-2f 



( I a{r, x) a(0, 0) | )| ~ K a )p ^ ^ e'^^^^) "\ (40) 

h 

It is readily verified that the weak field matrix elements p9|) lead to the 
dynamical structure factor: 



5 ( a; — 1 /m? + A;2 ) 
5(u;, k)=^h\{a)fY. — (41) 

b 

We illustrate this dependence in Fig. [51 The two-particle continuum of Fig. ^ 
has been replaced by a series of sharp peaks — the elementary excitations 
are now confined. 

In the next section we shall continue our discussion of the bound state 
matrix elements (jHEJ and ()39|) . We shall demonstrate how they (and indeed 
those in other models) follow from a direct computation within the so-called 
'two-quark' approximation [22]. In the same spirit as the relativistic mass 
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Figure 2. Ising model dynamical structure factor in the ordered regime and in a weak magnetic 
field. For illustration we have replaced the Dirac delta functions by unity. B = ■Kh\{cT)\'^ /m and 
we have set h = 0.001. 



corrections have been studied by Fonseca and Zamolodchikov [22] matrix 
element corrections may in principle be obtained. Their effect on the long 
distance asymptotics is presumably less significant however, and leading or- 
der results may be sufficient for many purposes. 



2.3. Bound State Amplitudes 

As we see from (|32() the amplitude we seek is determined by the overlap 
between the vacuum, the operator, and the bound state in its rest frame: 

(0|O|^f,). (42) 

The main problem is that we do not have an exact handle on the bound 
state l^'fe). However, following Fonseca and Zamolodchikov [22] for weak 
confinement (and £h = mfe — 2m <C m) we may consider the 'two-quark' 
approximation: 

1*,) = ^ r |^^,(<7)Zt(g)Zt(-g)|0), (43) 

where Z'^{q) are the Faddeev-Zamolodchikov creation operators of the con- 
stituent 'quarks'. In the case at hand these are fermionic creation operators, 
but in general they furnish rapidity dependent commutation relations. They 
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may also carry isotopic indices. The function ^b{Q) is the Fourier transform 
of the normahzed wavefunction (|27j) : 

/oo POO 
dxe*''^^b(x), where / fix = 1. (44) 

-oo J — oo 

In our conventions, the overall factor of l/y/m is required on dimensional 
grounds. An illuminating application of such a wavefunction to bound state 
mass spectra may be found in the work Gat and Rosenstein [25]; the inte- 
grable 1+1 dimensional massive Thirring model is viewed as a problem with 
delta-function confinement and the bound state meson is nothing but the 
sine-Gordon particle. In the present context, the approximation ()43() has 
the virtue that it embodies both the integrable aspects of the unperturbed 
system together with the most important effect of (non-integrable) confine- 
ment. As follows from the 'two-quark' approximation, the matrix element 
of the spin field is given by: 

(0|ct(0,0)|vI/,) = _ / dx^x) Jie^'i-{0\a{0,0)\q,-q), (45) 

where the normalized antisymmetric wavefunctions take the form: 

^bix) = ^ > 0. (46) 

' V2(mA)V6|Ai'(-Z,)r ^ ' 

Substituting q = m sinh 9 into the Ising form factor (|1H) one obtains 

( I a(0, 0)\q,-q)=i{a ) ^ (47) 

The Schrodinger description is valid for small momenta [22] and in view of 
this we expand the form factor in powers of momentum: 

(0|a(0,0)|g,-g)=i(a)A + o(^)'. (48) 



To lowest order equation becomes: 

{a) d^'b(x) 



(0|a(0,0)|^,)- 3 

1712 



(49) 



x=0 



Substituting in the explicit form of the wavefunction (|4f)|) yields the bound 
state matrix element: 

(0|cT(0,0)|*,) = (-l)^i^^/^, (50) 

m \ 2 
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in agreement with equation (|38|) . There are a couple of things to note in 
this relatively simple example. First, to lowest order the bound state de- 
pendence has disappeared: the normalization of the Airy wavefunction (|46|) 
has canceled under the differential action of the form factor H48() . Second, 
the derivative of a symmetric wavefunction must vanish at the origin, and 
by equation (|49|) would yield a zero matrix element. 

We note that whilst O {£b I'm) corrections to H50() are easily obtained from 
the expansion ()48|1 , an accurate determination of the coefficients presumably 
requires the Bethe-Salpeter approach of Fonseca and Zamolodchikov [22]. 
This lies beyond the scope of this current work. This therefore completes 
our brief discussion of the Ising model in a magnetic field. We turn our 
attention to another model of condensed matter physics which also exhibits 
aspects of confinement. 

3. Quasi-ID Wigner Crystals 

In this section we describe a possible realization of confinement in Wigner 
crystals [26]. If the average potential energy of an electron system exceeds 
its kinetic energy, the electrons themselves may adopt a regular arrangement 
in space and thus form a crystal. Wigner (or electron) crystallization is a 
spectacular consequence of electron-electron interactions. 

One may gain some insight into Wigner crystallization from a simple 
one-dimensional example — the extended Hubbard model at quarter-filling: 



c\ ^ is an electron creation operator and nj^o- = c[ ^Ci^a is the number operator. 
The electrons move on a rigid'^ lattice with sites labeled by i. They carry spin 
a =t, i and obey the usual anticommutation relations. The most important 
effects of the long-range ID Coulomb interaction are described by the on-site 
repulsion U and the nearest neighbor "extended" repulsion V|| . The model 
is defined to have on average one electron for every two sites; in view of the 
Pauli principle, the band is quarter-filled and kp = 7r/4. For related works 
see for example [27-33] 

For small values of U and V|j the system (|51() is a metal. Indeed, at 
quarter-filling, elementary band theory predicts metallic behavior. How- 





(51) 



"^For simplicity we neglect the possibility of lattice distortions or phonons. 
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ever, as the interactions are increased, this system exhibits a metal-insulator 
transition [27,28,34]. This T = quantum phase transition [7] is driven by 
8kp Umklapp processes. Such processes become relevant only in the pres- 
ence of sufficiently strong repulsion. The insulating state, which arises due to 
interactions rather than band filling, is an example of a Mott insulator [35] . 

A glimpse into the nature of the Mott transition and the associated 
Wigner crystal, is possible in the limit [/ — > cx) [34]. In so far as the 
charge degrees of freedom are concerned, the model with V|| = may be 
mapped onto a model of free spinless fermions. Including the nearest neigh- 
bor repulsion one obtains a one-dimensional model of interacting spinless 
fermions [34]: 

He = -t ^(c1q+i + h.c.) + V|, ^ im - 1/2) - 1/2) , (52) 

i i 

where ui = "^^rii^^. As is well known [4,7] this model maps on to the 
spin- 1/2 XXZ spin chain via the Jordan- Wigner transformation: 

Hf'''^ = jT.(SfSf+^ + + A 5f 5f+i). (53) 

i 

The parameters are related by 

J = -2t, A = -^. (54) 

With repulsive fermions both J and A are negative. With this sign of J, it 
is known that for — 1 < A < the chain ()53() is gapless, whereas for A < — 1 
it develops both a gap and long-range antiferromagnetic Ising order: 

{{-lysny^O. (55) 

The gap implies that the spinless fermion model 1)52(1 is a Mott insulator 
for V|| > 2t, and a metal for V|| < 2t. That is to say, as U ^ oo in our 
original model (|51|) the Mott transition occurs for Vj| = 2t [34]. This is 
confirmed by numerical simulation [28] . As follows from the mapping to the 
XXZ chain, this metal-insulator transition is of the Berezinskii-Kosterlitz- 
Thouless (BKT) type. It also follows that close to the transition (when 
the insulating gap is much smaller than the bandwidth) the model is SU(2) 
invariant. At low energies we may also "linearize the fermionic spectrum" 
[4, 36] and 1 -|- 1-dimensional Lorentz invariance emerges. Away from the 
transition the SU(2) symmetry is broken down to U(l) x Z2, and the Lorentz 
invariant limit is no longer justified. We note that as U is lowered, one 
expects this critical value of V\\ to increase. Indeed, the numerical simulations 
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indicate that the phase boundary connects the "Ovchinikov point" {U, = 
(oo,2t) and the "Penc-Mila-Zotos point" {U,V\\) = (4t,oo) [27,28] — see 
Fig. 131 It is plausible that model (|52|) may also describe the transition for 
finite U > 4t, providing V\\ is replaced by an effective repulsion V^^^{U). 
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Figure 3. Metal-Insulator transition in the ID extended Hubbard model and the associated Mott 
insulating Wigner crystal or Akp CDW state. The phase boundary connects the "Ovchinikov 
point" (O) to the "Penc-Mila-Zotos point" (PMZ). Along the "XXZ line" (U — » oo) the system 
undergoes a BKT transition at the point O . 



Since Sf 



rii 



1/2, the ordering (|55() of the XXZ chain implies that the 



density of spinless fermions 



2)) alternates from site to site (see Fig. |3]) and 

((-l)Vi)/0. (56) 



i.e. the insulator is a Wigner crystal or ^^Akp CDW" (Charge Density Wave). 
Indeed, the interaction term in (|52() clearly favors this alternation. As shown 
in Fig. 13 the ground state is doubly degenerate with respect to charge or- 
dering. The elementary excitations or 'quarks' in this insulating regime are 
domain walls separating these two ground states. It may be seen that these 
'quarks' carry fractional charge ibe/2. Moreover, it may be shown that these 
'quarks' are deconfined. As in the Ising model, this manifests itself in the 
absence of poles in the appropriate dynamical (charge) response functions. 
In the case at hand these are the spectral function [31, 32] or the optical 
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Figure 4. Degenerate Wigncr crystal (or CDW) ground states in tlic Mott insulating regime. 
The darker regions are sites of higher electron density. The elementary excitations, or 'quarks', 
are domain walls separating these ground states. 

conductivity. The relevant probes couple to pairs of 'quarks'. 

This concludes our discussion of the purely one-dimensional system (|51|) . 
In the next section we shall discuss what happens when weak interchain 
interactions are switched on. 

3.1. Interchain Interactions 

Let us take our model of spinless fermions (|52)) and switch on weak interchain 
interactions, VI ^ Vj| : 

He = -t^{4j Ci+i,i + h.c.) + I^i Y.i'^iJ - V2)(ni+i,j - 1/2) 

+V± 5^K, - l/2)K,+i - 1/2). (57) 

ij 

For notational simplicity we consider a two-dimensional array of chains. The 
extension to three-dimensions is straightforward. We also neglect the possi- 
bility of interchain hopping. In the ordered phase, we may treat the inter- 
chain interactions in a mean field approximation. That is to say, from the 
perspective of a single chain, we replace the neighboring chains by a suit- 
able expectation value. This yields a set of decoupled chains in a staggered 
chemical potential: 

He = -t ^(cJq+i + c\_^iCi) 

i 

+V\\ i^i - 1/2) (ni+i - 1/2) + 5± J^(-l)' n,, (58) 

i i 

where 



(59) 
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and is the transverse coordination number. As follows from H56|) the 
expectation value occurring in (|59|) is non-zero even for Vj_ ^ 0. This 
interaction renders the metal-insulator transition first order. 

We shall discuss the model (|58() in the limit where the dynamically gen- 
erated 'quark' mass is much smaller than the bandwidth, m <^ t. In this 
limit, close to the transition, Lorentz invariance emerges and field theory 
may be applied. In the vicinity of the SU(2) invariant Ovchinikov point 
the mean field Hamiltonian H58|l may be bosonized [36]. The corresponding 
mean field" action may be written: 



1 - /3 ' 

— {df,<^cf - A cos /3 ^>c + A cos - $c 



(60) 



where in this normalization = 1. This model is the well studied one- 
dimensional double sine- Gordon model — see for example [37,38] and refer- 
ences therein. The coupling constants ft and A are related to the parameters, 
U,V\\,V±,t, of the bare Hamiltonian. In particular, A^ ~ (^|/2i — !)• The 
term ft cos /3 $c represents the 8A;i?-Umklapp processes, and A cos ^ $c repre- 
sents the interchain interactions. We note that the latter term may also arise 
from purely ID intrachain dimerization [31]. Indeed, the confining aspects of 
the double sine-Gordon model were first applied in Affleck's insightful study 
of the Heisenberg chain with frustration and dimerization [2]; for related 
works see also [39-41]. 

The model H60() . with = 1, has two characteristic soliton mass scales 
generated by the intrachain and interchain perturbations: 



A^~t/i'/'e-i/A, (61) 
M = A3^ = Cp/^ (62) 

where C is given by [42,43]: 

2 r[i/6] /^r[3/4]y/^ 

V^r[2/3] V2r[i/4]; • ^''^ 



The ratio of these scales, M/m, may be taken as a measure of the strength 
of the confinement. In particular, as we shall discuss below, the number of 
stable bound states decreases as M /m is increased. 

Perhaps the simplest regime to consider, is that with M ^ m. In this 
limit, m and fi ^ and (|60|) yields the = 1/4" sine-Gordon model. The 



"The most essential corrections to the mean field action couple neighboring chains n and m via 
an interaction ~ V± J dx^'^^^^c^ ■ This lifts the SU(2) degeneracy and also loads to an 
anisotropy in the optical conductivity tensor cr^^ ^ ct" ^ . 
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spectrum of this model is known to consist of a triplet of mass M and a 
singlet of mass \/3M [2,3,43,44]. As such, the spectrum in the charge sector 
consists of four bound states. The three lightest 'mesons' are degenerate 
in mass, and are distinguished by their charges, ite, 0. In the language of 
the XXZ chain, these charges correspond to the projections of the auxiliary 
spin variable 5^; we emphasize that these particles are actually spinless. 
Their properties are summarized in Table ^ and in Fig. |S1 The relatively 



Table 1. Bound states 



Mass, M 


Charge, e 


Till = 1 


+1, 0, -1 


m2 = \/3 






small number of bound states in this regime parallels the reduction to the Eg 
spectrum in the h ^ Ising field theory as m — > [20] . The bound states will 
manifest themselves as sharp peaks in the dynamical response functions. In 
particular, the optical conductivity (as determined by the current operator) 
will probe the neutral bound states. Their contributions may be extracted 
using the known sine-Gordon form factors. 

The regime M <C m bears a close relationship to the Ising model in a weak 
magnetic field. In this regime, the confining perturbation is very weak and 
one may (in principle) apply the Bethe-Salpeter approach of Fonseca and 
Zamolodchikov [22]. Although we have not performed this rather technical 
analysis for the double sine-Gordon model, the lowest order Schrodinger 
description is applicable: 

^ +\\x\^ix)=£b^bix). (64) 



m 



Solitons (whose size is much smaller than their separation) are treated as 
point-like particles interspersed by a region of false vacuum. This gives rise 
to a linear potential or 'string'. The string tension A is related to the bare 
tension A via: 

A = 2A (cos(/3$c/2)) . (65) 

It may be seen that A has dimensions of [mass]^ and (|64() is consistent on 
dimensional grounds. Note that quantum mechanical fluctuations are incor- 
porated in the expectation value appearing in This quantity behaves as 
and reduces A with respect to its classical value [2] obtained as 0^ 0. 
As in the weak field Ising model, the spectrum consists of a plethora of bound 
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states, with masses given by: 



TTift = 2m + 



m 



2 \ 1/3 



+ Ai(-Zb) = 0. 



(66) 



This formula apphes to both charged and neutral bound states, and the 
number of distinct masses is determined by the stability threshold, < 
2m\. In this regime we therefore have many generations of 'mesons' indexed 
by the zeroes of the Airy function — see Fig. [HI Each generation is four- fold 
degenerate and consists of a triplet with charges ite, 0, and a singlet with 
charge 0. This degeneracy is evidently lifted as the ratio M/m increases. 
The evolution between these limits is depicted qualitatively in Fig. |S1 The 
contributions to dynamical response functions may be calculated along the 
lines of Q 



a 
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Figure 5. Qualitative evolution of the excitation spectrum in the weakly coupled Wigner crystal. 



Having discussed the effects of confinement, let us briefly digress on the 
effects of a finite on-site Coulomb interaction on the spin sector. At finite 
U there are two effects. First, there will be a Heisenberg exchange between 
spins on the same chain. Second, there will be a spin dependent contri- 
bution to the interchain interactions arising from the Coulomb repulsion; 
in the bosonization approach, the electron density operator contains a 2A:f 
component which is spin dependent — see for example chapter 30 of ref. [4]. 
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In this way one may derive an interaction of the form: 

F""^ = J cos («>^/2) cos ($™/2) cos ((^^ - /4) . (67) 

Since cos(<&s/2) is the bosonized version of the staggered energy density for 
the Heisenberg chain, and spins at 1/4-fining are next-nearest neighbors, it 
is plausible that 

y"*" = J' (Si.S.+2)n(S..Si+2)™ cos (($^ - $^)/4) (68) 

is also valid when the spin exchange is much smaller than the charge band- 
width. This interaction leads to spin-Peierls ordering at low temperatures. 



4. Conclusions 

In this work we have discussed two instructive examples of confinement 
which occur in low-dimensions. Although such problems are generically non- 
integrable, field theory approaches are able to shed light on many interesting 
aspects. In closing, let us touch upon a somewhat more exotic form of con- 
finement in a model with two scalar fields: 

i=l * 

where 

Vc = -92 cos $2 cos(^>i/4). (70) 

First, note that $2 is massless in the absence of the interaction Vc- There- 
fore, unlike the previous examples considered, the expectation value of the 
confining operator H7U|) vanishes in the unperturbed model. Since the mass 
scales for the <l>i and <I>2 excitations are given by 

1/2(1-/3?) / /32/8\ V2(l-/3|) 



the corresponding string tension A ~ 771-2 ^ nonlinear dependence on 
the coupling g2. Second, let us consider a pair of $1 solitons interpolating 
between (0, 27r) and (27r, 47r) respectively. In this background, the operator 
cos($i/4) changes sign, and the interaction H7U|) promotes massless fluctu- 
ations of ^>2 in the region between the solitons. In this way it is possible to 
argue that the potential between the solitons is of the form 

U{x) = X\x\ + (72) 
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where Z{nk) is a function of the occupation numbers of the gapless modes. 
In this case, the $i solitons are confined by a double-well potential. 
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